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ABSTRACT 



The relation of the total structure weight of a 
two-dimensional ring frame to the stress distribution in each of 
the elements is investigated. The total weight of the two- 
dimensional structures investigated were a minimum when design 
stresses had been attained in all elements of the structure. 

An expression is developed that approximates the change 
in stress distribution in a member when additions or deletions of 
material are made in the member itself and its associated members. 

This expression is used as the basis for a rational 
design technique that converges to a ring frame of minimum total 
weight. 
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I. INTRODUCTION 



The transverse structure of a wall sided ship consists 
primarily of bo t tom and deck beams, side frames and stanchions. 

This structure is called the transverse ring frame. The loads this 
structure is subjected to are side and bottom water pressure and the 
ship weights as well as dynamic loadings. 

The structural analysis of this type of structure, in two 
and three dimensions, is discussed in considerable detail in the 
literature. Computer solutions for both the two and three dimensional 
aspects are available. These rapid analysis tools are necessary pre- 
requisites for a rational design synthesis. 

In its most general form the optimization of the transverse 
structure entails the specification of the transverse frame spacing, 
the effect of longitudinal frame spacing, and the distribution of 
material throughout the frame. This thesis deals with the last part 
of the process, the optimium distribution of material throughout the 
ring frame of a wall sided ship once the transverse and longitudinal 
frame spacings are fixed. The optimum structure is defined as the 
least weight combination of deck beam, side frames, etc,, that will 
carry the design loads. 



- 1 - 



II. MINIMUM WEIGHT STRUCTURES 



The end and field moments that a loaded beam element develops 
are directly related to its end restraint. These moment distributions 
are subject to wide variations. The bending moment diagrams for a 
beam element supporting a distributed load are reproduced in Figure 1, 
Full rotational restraint at both ends produces a symmetri- 
cal moment distribution and one in which the maximum moments occur 
at the ends. The field moment dominates when the moment restraint 
at the ends is zero. The intermediate restraint case is character- 
ized by a symmetrical moment distribution and identical end and 
field moments. In every case, limiting or design stresses are 
reached only at discrete points in the element. 

The size or weight of any element depends on the maximum 
moment developed when the design load is applied. The wide variation 
in maximum moments is accompanied by corresponding variations in the 
weight per unit length of members. Table 1 lists, for the ideal 
section developed in Appendix C, the variation in weight for the 
structural configurations of Figure 1, The absolute minimum weight 
per unit length occurs when the maximum stress occurs at three points 
simultaneously in the beam element, The specific weight increases 
by thirty percent going from the fixed-fixed beam to the free-free 
or fixed-free beam. Doubling the maximum moment for this ideal 
section increases the section weight by sixty percent. 
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The weight ratios of Table 1 underline the fact that the 
total weight of any structural system is related to how the individ- 
ual elements interact, A minimum total weight of structure can be 
obtained by substituting members whose moduli satisfy the inter- 
mediate restraint condition. In most cases, the actual moment 
distributions that obtain when the structure is assembled will not 
be identical to those assumed conditions. The final total structural 
weight will depend on how the material distribution has been altered 
to obtain convergence. This, in turn, depends on the design criteria. 
The examples that follow illustrate how the total weight varies when 
the design criteria is that design stress must exist in at least one 
element; the stress in all other members must be equal to or less than 
the design stress. 

The first example. Figure 2, contains only elements sub- 
jected to bending. The length to depth ratio cf each element is 
assumed to be large enough to exclude shear effects. 

An initial choice of members based upon the intermediate 
restraint condition and whether the elements carry an external load 
gives : 

1/16 W hi 

h - ; z i = 0 = z 3 

Design 

The structure made up of these members does not possess adequate 
strength. The maximum stress in member two is greater than the 
design stress. The minimum weight members do not form an acceptable 
structure. 
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This structure cAn only be brought up to an acceptable 
level of strength by adding material, The weight of the final 
structure satisfying the design criteria will vary considerably de- 
pending on where the material is added. There are two directions 
available. One is to add material so that the moments at the nodes 
are very small. This will cause the field moment in member two to 
predominate, The other members will be very small. The second 
approach is to increase the size of all the members. This will 
bring about redundant moments at the nodes which are substantial, 

The weight curves of Figure 3 illustrate the trade off 
possibilities that exist when all the members are varied. Large 
increases in the specific weight of members one and three are 
accompanied by relatively small increases in member one as the dis- 
tribution factor approaches one. These effects are summarized in 
the total weight curve, These curves , for all the cases looked at* 
have definite rrinimums. When “ 5L^ = 5L^, the curve shown in 
Figure 3, the minimum total weight occurred when all member sizes 
were equal and the node moments were maximum. At this minimum point, 
design stresses were obtained in all members. 

The following example is identical to the preceding except 
that all members have to support an external load. The loading is 
such that the fixed end moment of the middle span is twice that of 
the adjacent members; load symmetry is maintained. 
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As before, the structure made up of least weight members, 
corresponding to the intermediate restraint condition, is not 
satisfactory. All members are cverstressed , Sufficient strength 
can only be obtained by adding material. 

Total weight curves for this structure are shown in 
Figure 4, Each of these curves have minimum points which occur at 
k * .5. At this point, all members are of equal size and the node 
moments govern. It is obvious then that the minimum total weight 
occurs when design stress levels exist in all the members. At all 
other values for the distribution factor, one or two members have 
too much material. 

Only bending stresses were considered in these examples. 
The substitution of a vertical member at each of the supports re- 
quires that direct stresses be accounted for (see Appendix B) * The 
total weight curve for this structure also has a unique minimum. 
This minimum occurs when design stresses exist in all members. 

The iteration technique developed in this thesis conver- 
ges on a structure which, when the design load is applied, develops 
design stresses in each member simultaneously. These stresses may 
develop at the ends of the members or at any point along the 
members. The structures which this design process converges to 
represents the optimum distribution of material for the given loads 
and the fixed frame spacing. The total weight of these structures 
is a minimum. Any other combination of members will result in a 
higher total weight of structure; the stress in some or all the 
members will be less than the design stress. 
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MOMENT DIAGRAMS FOR A BEAM SUPPORTING A DISTRIBUTED LOAD 



FIGURE 1 
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TABLE 1 



Maximum Mome nt 
1/16 \TL 7 
1/12 WL 2 
1/8 WL 2 



Posit ion 
End, Field 
End 

End or Field 



Mo neat Ratio 



M 

o 

1.33 M 

o 

2.0 M 

o 



Weight : Ratio 



W 

o 

1.21 W 

o 

1.59 W 

o 



WEIGHT AND MOMENT RELATIONS FOR THE 
BEAM SUPPORTING A DISTRIBUTED LOAD 
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STRUCTURE FOR EXAMPLE 1, SECTION II 




STRUCTURE FOR EXAMPLE 2, SECTION II 



BEAM SYSTEMS FOR EXAMPLES ONE AND TWO 



FIGURE 2 
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Total 

Structure 

Weight 
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Weight Per 
Unit Length 



Total Weight Curve 




Distribution Factor 



TFtE TOTAL WEIGHT AND WEIGHT TER UNIT LENGTH 
CURVES FOR EXAMPLE ONE, SECTION II 



FIGURE 3 
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Total 

Weight 




.5 



1,0 



Distribution Factor 



THE TOTAL WEIGHT CURVES FOR EXAMPLE TWO, SECTION II 



FIGURE 4 
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III. CONVERGENCE PROCEDURES 



The properties of a two dimensional transverse ring 
frame of a wall sided ship ray be considered similar to those of a 
plane frame. The members of this plane frame are subjected to bend- 
ing and direct stresses; side sway is not considered. In its most 
complex form, this space frame is made up of many bays and decks. 
Picking member sizes that are adequate, i.e, possess 
enough strength to support the applied loads, can be accomplished in 
a straight forward manner, St. Denis [1], presents a rational 
approach to this problem. The refinement of these initial* choi ces , 
i.e. cutting the strength of each member until it just meets the 
design stress, is discussed in this section, 

A two dimensional space frame of more than one deck is a 
statically indeterninent structure. The stress distribution result- 
ing from the design loads depend upon the size of each member of the 
frame. Changes in these stress distributions are expected if the 
size of one member is altered. The effects of changing the size of 
a member in a reasonably balanced structure will be greatest in the 
immediate vicinity of the member. It is assumed in the following 
that the stress or moment distribution in a member is affected only 
by changes in members that join the defining nodes of the member. 
This is the same as saying that each member of the frame is the 
same as member 2 of Figure 5, The stress distribution in member 2 
is affected by size changes in any of the members shown. 
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The refinement problem is this: given a frame with a fixed 



layout and subjected to specified design loads, what distribution 
of material will produce stress distributions such that in each 
member the maximum stress equals the design stress. This means 
that in each member that has a maximum stress less than design stress, 
the section is too heavy. Material must be taken away or the stress 
distribution altered so that the member carries more of the load. 

If the maximum stress in the member is greater than design stress, 
the member is undersized and the opposite holds. 

An approach to the problem that is suitable to machine 
applications is a selective procedure. After each stress analysis, 
the member in the worst condition, over or under design stress, 
is selected. The size of the member is altered based upon existing 
maximum moment. The structure with the corrected member undergoes 
another stress analyses. This cycle is repeated until all members 
are at the design stress level, 

A more appealing approach is one in which changes can be 
made in all the members before each stress analysis. The results 
of each stress analysis are used as the basis for changing unsatis- 
factory members. 

The implementation of this series approach requires that 
interaction effects be estimated. Unless account is taken of how 
the stress distribution changes when the distribution of material is 
altered, the system will be subjected to undesirable oscillations. 
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The member shown in Figure 6 depicts a general member taken 
from the two dimensional plane frame* It carries a distributed 
load, W, and has non-syrnmetrical end restraint* The stress dis- 
tribution is not symmetrical; node one is assured to be overstressed. 

The information required is how will this moment distribution 
change when any or all of the members involved change* At each, node 
three conditions have to be considered, an increase, decrease, or no 
change in the joint rigidity. 

An increase in the size of all the members adjacent to node 
two will provide additional rotational restraint. This is accompan- 
ied by an increase in the moment at node two and a decrease in the 
maximum moment. If the changes in the members adjacent to node two 
were sufficiently large, no change in member one would be required 
due to the altered stress distribution. 

Increasing the size of member one alone will cause the 
maximum moment to increase. If the size of the new member is picked 
based on the existing maximum moment, then it also will be over- 
stressed when the structure is reassembled. 

Table 2 is a summary of the conditions that have to be 
accounted for. Column three has the information needed to make a 
rational choice of new members. An expected increase in the maximum 
moment indicates that a larger than usual correction should be made 
in the already overstressed member. This would be accomplished by 
adding an increment to the existing moment," The new member size would 
reflect this added increment. 
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Additional tables would be required for each of the 
following: a field moment governing, end moment identical and govern- 
ing, and non loaded elements. 

The Hardy Cross method of moment distribution offers a 
direct method for obtaining numerical estimates of changes in maxi- 
mum moments due to shifts in material distribution about nodes. The 
expression that is derived below actually embodies the tables dis- 
cussed in the last section, 

A first approximation to the moment U of Figure 5 can 
be. obtained by carrying out the Hardy Cross method of moment dis- 
tribution through the first unlocking of all eight joints. It is 
not necessary to require that the joints three through eight be 
fixed , 

First release the joints adjacent to node two and 
distribute the unbalanced moments, U, The carry over from joints 
four, six and eight to node two is: 

CO - - 1/2 k 62 U 6 - 1/2 k g2 U 8 - 1/2 k 42 U 4 

This carry over is added to the unbalanced moment at node two, . 

The amount of this summation that carries over to node one after 
node two is unlocked is: 

1/2 k 21 [1/2 (k 62 U 6 + kg 2 U g + k A2 U^) - U 2 ] 

The carry over from the joints adjacent to node one is 
found in the same way. It is: 

CO - - 1/2 k 51 U 5 - 1/2 k 71 U ? - 1/2 kg 1 Ug 
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'The moment of node one in member 1-2, M rf% can be 
expressed in the following approximate form: 



M 
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F[ * [ l-2 " k 12 C 1 + 111 k 12 ^ k 51 U 5 + k 31 U 3 + k 71 U 7 > 



1/2 U 2 + 1/2 k 12 k n U 2 + 1/2 k 21 (k 62 U 6 + 



k 82 U 8 + k 42 V 



1/4 k 21 k 12 (k 62 U 6 + k 82 U 8 + k 42 U 4 > + 1/4 k u k 21 U 2 



1/4 k i2 k 21 U 2 



( 1 ) 

If M is the maximum moment in the member, the change in 
the moment as the adjoining members and member two itself changes is 
desired. The quantities FEM^ and depend only on the loads and 
the geometry and therefore are constants. The expression for 
can be looked upon as a function of the distribution factors, , 
which in turn depend directly on the section properties. 



M 12 ” f (k 12 k ij ) 

The change in is 

tw 3f . 3f 

dM i2 - "ak~ dk i2 + ie: 

12 ij 

AM = av + — ^ — 

12 3k. _ 12 3k., 

12 ij 



or 



dk. . 
ij 



Ak. . 
ij 
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Applying this to equation (1) results in: 



AM 



12 



- (q ik 12 + 1/2 (k 51 U 5 + k 31 U 3 + k n U ? ) «, 12 



+ 1/2 k 12 <ik 5X U 5 + 4k 31 U 3 + 4k 71 V - 1/2 U 2 4k 21 



+ 1/2 U 2 (k 12 Ak 21 + Ak 12 t' 21 ) + 1/4 (k 62 U & + k g2 U g + 



k 42 V Ak 21 



+ 1/4 k 21 (Ak 62 U 6 + Ak 82 U 8 + Ak 42 V + 1/4 U 4 



( k i2 Ak 21 + Ak 12 k 21 ) 



( 2 ) 

An adequate comparison of the orders of magnitude of the 
different parts of this expression requires knowledge of the joint 
unbalance. If it is assumed, as it was previously, that the loads 
are reasonably balanced throughout the structure, the structure of 
Figure 5 is again applicable. Applied to this structure equation 
(2) reduces to: 

AM 12 * ~ U i Ak i2 ~ 1/2 U 2 Ak 21 + 1/2 U 2 (k 12 Ak 21 + Ak 12 k 21 } 

+ 1/4 iq (k 12 Ak 21 + Ak 12 k 21 > 

(3) 
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This is the same as neglecting the unbalance at joints 
three through eight. The inclusion of these terras in a computer 
program is straightforward, but of course, would require more 
computational time. This final approxinateion accounts for changes 
in any of the seven members involved, A more convenient form for 
equation (3) is: 



AM. . * [-U. (1 - 1/4 k .) + 1/2 k. . U . ] Ak. . 
1J i ji Ji J ij 



+ 1/2 [- U. (1 - k^.) + 1/2 U. k..] Ak.. 

J ij i ji 



(4) 



Equation (4) gives a numerical value for the expected 
change in M if the adjoining members and member itself change size. 
This is an approximation as to how the moment distribution in member 
ij is affected when any or all the members change size. This is a 
valid approximation for both loaded and unloaded elements. There 
is no restriction to members undergoing only bending stress; it is 
equally useful for members with combination stresses. If different 
types of materials are used, their effects are included in the 
distribu tion factors , 

The use of this expression to find changes in the field 
moment is a little more involved, Equation (4) has to be evaluated 
for each end of the member in question if the moment distribution 
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is not symmetrical and added to the existing end moments, Ihese 
new moments combined with the existing end shears give an approxi- 
mation to the change in field moment, 

A word about signs is appropriate. Clockwise moments 
are positive , counterclockwise moments are taken to be negative. 

The signs that result when AM„ is calculated cannot be interpreted 
literally, A positive value for AM_ may indicate an expected de-‘ 
crease in the moment, * It is necessary to check the sign of 
the existing moment and then interpret the sign of AM , 
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A SUBSECTION OF THE TW0-DI2IEIISI0NAL RING FRAME 



FIGURE 5 
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A GENERAL BEAM MEMBER OF THE TWO-DIMENSIONAL RING FRAME 



FIGURE 6 
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TABLE 2 



Cumulative Effects of Member Variations 



End 

w/ 

Mom) 



Min 



End 

w/* 

Mom) 



Max 



( 1 ) 

Adjacent 

Members 



+ 

o ‘ 



+ 

o 



( 2 ) 

Member (3) 

to be Exp, Change 

Changed in Max. Mom, 

+ Increase 



4- Nc Change 

+ Increase 



4* Decrease 

4* No Change 

4- Decrease 



Explanation of Table 

This table is drawn up for Member Two of Figure Five. It is 
assumed that the moment distribution is nonsynmetrical and that 
Member Two is overs tressed , Column three gives the expected change 
in the maximum moment based upon the cumulative effects of member 
changes summarized in columns one and two. 

Column 1: The sign indicates how the overall stiffness of 

the joint is expected to change based upon the 
stress conditions in the adjacent members. 

Column 2: The size of Member Two is expected to increase 

Column 3: The expected change in the maximum moment based 
upon the data of columns one and two. 
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IV. THE DESIGN SPIRAL 



The result of the previous section is a technique for 

estimating moment changes. These moment changes are brought about 

when the members at a node change their size. 

This is of use in the design process because it presents 

a method of estimating not only the magnitude of changes in the 

stress distribution (maximum moment), but also the direction of 

the changes, It permits a more rapid convergence to the final 

structure because each member of the structure can be modified (if 

>*■ 

necessary) on each pass. And also because the members are picked 
based upon conditions that will prevail; not those of the last 
stress analysis. It is essentially a method of predicting the 
final moment distribution and picking the members based on this 
distribution. If the moment changed predicted were exact instead 
of approximations there wouldn ? t be any iterative procedure 
necessary, 

A desirable aspect of this technique is that it is self 
modifying when applied repeatedly, The predicted moment changes 
are updated each time a new member is picked. In this way oscilla- 
tions are avoided and smooth convergence is assured, It is not 
necessary to perform a new structural .analysis until this self 
modifying aspect is completed. This occurs when there are no 
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significant* differences in the predicted moment changes from one 
cycle to the next. 

The rational design process developed in the sections that 
follow converges to the optimum distribution of material for a trans- 
verse ring frame. The general procedure followed is to use the 
results of a stress analysis as the basis for modifying the material 
distribution in such a manner that design stresses will be attained 
in each member. The results of Section II indicate that a structure 
which develops design stresses in each member simultaneously is 
the least weight structure, 

>■ 

Equation (4) is used in this process to make corrections 
to the results of the stress analysis, A change in the scantlings 
of any member of the ring frame changes the stress distributions. 

As pointed out in Section III and summarized in Table II, the 
effect of a change in scantlings depends not only on the magnitude 
of the change! but also the existing moment distribution. Once, the 
maximum moment in the member is identified, equation (4) gives 
direct information as to how this moment will change when the 
scantlings of the members of the ring frame are revised. 

This correction to the existing moment distributions 
provides the means of facilitating rapid convergence to the optimum 
structure. Changes in member scantlings are based on these pre- 
dicted moments. The objective of design stress in each member of 
the ring frame is obtained much more readily. 
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■ The use of equation (4) in the design process requires 
the standard output of a stress analysis* end forces* plus end 
moments. The process begins with the results of the stress analysis 
which used the initial choices of members. 

The computer begins with member one. The governing 
moment is ascertained and the stress level is computed based on this 
maximum moment. If the stress level is within a certain prescribed 
range of the design stress* the member is not changed. This pro- 
vides a means for ending the iterations. If the stress level is 
not satisfactory, a new member size is calculated and stored in the 
matrix* K, This procedure is repeated for each member of ''the 
structure. The result of this cycle is the matrix K. The members 
that make up K represent the first cut at bringing about the desired 
distribution of material. These members are fictitious in the 
sense that they are not used in an actual stress analysis. The 
properties of these members are based upon the existing stress 
distributions. 

The computer now’ starts to revise the members of matrix K; 
it starts with the first loaded member. The stiffnesses, k..,of 

ij 

the members at the defining nodes of the member to be changed are 

computed using the data of matrix K, This data is used to compute 

AM,, or AM., or both depending on the stress distribution. If the 
ij Ji 

field moment governs AM^ is computed for each node; when an end 
moment governs A!I_ is calculated for the governing node. 
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The maximum moment of the member to be changed is modified 
by adding the correction factor A>k_. . The properties of a new 
section are calculated based on this modified moment. These 
properties are inserted in the proper slot of matrix K , 

The computer chooses the next loaded member and repeats 
the above. The unleaded members are revised after new properties 
have been picked for all the loaded members. This requires fewer 
steps because the end moments always govern. 

The result after all the members have been looked at is 

» 

a new- matrix of members , K f A stress analysis could be initiated 
at this point. The results of the examples indicate, however, that 
a smoother convergence is obtained if this process is repeated, i.e, 

it 

find the matrix K , 

A stress analysis is now undertaken using the members of 

n 

the matrix K , The results of this stress analysis are used as 
the starting point for a new design cycle. 

Flow charts containing the details of the design spiral 
are contained in Appendix A. Illustrations of the design spiral 
in action are contained in Appendix B. 



V. RESULTS AND CONCLUSIONS 



A rational design procedure has been developed which 
permits the rapid synthesis of the transverse ring frame of a 
wall sided ship. The resulting transverse structure represents 
the optimum distribution of material, the minimum wieght combus- 
tion of deck beams, side frames, etc., within the constraints of 
fixed longitudinal and transverse frame spacing. 

Detailed applications of this procedure are carried out 
in Appendix L. Each step follows the subroutines developed for 

s- 

computer adaptation of this technique. Mo attempt has been made 
to adapt the procedure to hand calculations. 

This method converges rapidly. Large member changes were 
required in each of the examples to attain the design stress in each 
member. Example I required only one additional stress analysis 
after the initial choice of members were analyzed; Example II re- 
quired two additional moment distributions. 

The speed with which a satisfactory ring frame is obtained 
depends directly on the following: the limits set within which the 
member stress must fall and the use of weighting functions in con- 
junction with the predicted moment changes. 

The limits set for member stresses in Appendix B were 
liberal. This, of course, was conducive to a rapid convergence. 
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These limits are not fixed, however, and can be set at any desired 
level, It is possible to put new limits on the members stress with 
each problem. 

The predicted moment changes in the examples were used 
as given by equation (4) , The true approximation in equation (4) is 
in the magnitude of expected moment change. The direction of the 
change, plus or minus, is exact if the change is greater than a few 
percent. In all the cases investigated, the magnitude of the change 
was always low. Convergence would have been more rapid if the 
changes were multiplied by a weighting factor. This indicates that 
weighting factors, if necessary, can be developed as experience is 
gained with the basic program. These weighting factors would become 
part of the basic program. 

The results of Example II of Appendix B indicates that 
members whose continuity is not required to carry the applied loads 
may be eliminated. This is easily avoided by applying nominal 
loads to these members, e,g, their own weight. In an actual 
structure, this will seldom be a problem. 

The use of the ideal section of Appendix C was a necessity 
in the development stages of this program. The implementation of 
this design procedure, i,e, the actual programming, does not re- 
quire this ideal section. If a tape of actual sections is available, 
the program can use it and still work, The effect this would have 
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on the speed of convergence would be mixed. The actual members in 
most cases have larger section moduli than those the program would 
pick. The use of actual sections requires,ir, addition, the develop- 
ment of a trial and error solution in the subroutine COMSTRrSS , 
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APPENDIX A 



An algorithm covering the complete design process is 
developed, Flow charts covering the details of the individual steps 
of the algorithm are also presented. 



Design Algorithm 
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The flow charts are made up in the forn of subroutines. 
They cover all the steps presented in the design algorithm. 

The first flow chart is subroutine Moment, It uses as an 
input the data of the preceding stress analysis. This subroutine 
primarily finds the maximum moment and calculates the maximum 
stress. When field moments govern the member is given a special tag 
Subroutine Change checks the maximum stress in each member 
It stops the iteration process if all members are satisfactory. 
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Subroutine Stiffness uses the data of matrix K and the 



data of initial member choices to calculate the expected changes in 
stiffness. It also determines the moment change and adds this to 
the maximum moment, The nex; member parameters are obtained by 
calling the appropriate subroutine, i,e, Pick cr Comstress, This 
subroutine assumes that a member numbering system is available which 
facilitates picking the required data from the matrix K, etc. 

Subroutines Pick and Comstress both choose new member 
properties, Comstress handles members that have significant direct 
stress as well as bending stresses. 

The calling sequence for these subroutines x% T ould be: 

CALI. SUBROUTINE MOMENT 
CALL SUBROUTINE CHANGE 
CALL SUBROUTINE STIFFNESS 

CALL STRESS - or any suitable stress analysis program. 
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SUBROUTINE MOMENT 
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SUBROUTINE CHANGE 
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SUBROUTINE STIFFNESS 
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SUBROUTINE PICK 




SUBROUTINE CO'tSTFESS 
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appendix b 



mmEs 



EXAMPLE I 



The seven ner.be red frame shown on page 45 is used to 
illustrate the optimization technique developed in the body of 
the thesis, 

A, S e 1 ec tip n o Initial Member Proper tie 

' The initial set of member properties is part of the data 
input to the computer. The following were selected for this 
example. 

(a) Z 2 = 240 in 3 ; I = 2000 in 4 

(b) Z 1 = Z 3 = 120 ir 3 ; I = I = 820 in 4 

(c) Z ■ Z = Z = Z = 120 in 3 ; I = 820 in 4 

4 5 6 7 4 

These properties were obtained using the ideal section of Appendix 

B. < Stress _ Anal y_sj.s^ 

The Hardy Cross method of moment distribution was used to 
obtain end moments and reactions. The results of this first analy- 
sis are: 



M 



12 



- = M- ) . t = 684 KIP FT 



21 



Max 



M 17 M 28 " M 2 1-fax “ 495 KIP Ff 

M 13 = Vllax = 95 K ’ ; F 4 = 130 KIP FT 
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C, Subroutine Chance 



(1) The maximum stress in each member is compared to 
the design stress = AO KSI 

(2) The size of each unsatisfactory member is changed , 

The maximum stress in each member is; 

o-).. = 3A KSI 

1 Max 

a.,).. = 50 KSI 

7 Max 

a ) = 15 KSI + 10 KSI = 25 KSI - a combined stress member, 

3 Max 

The new properties for each member are: 



_ 2 Max . 3 T irnri . A 

7 = E 20 A m ; I„ ~ 1600 m 

2a. ’2 

Design 

Z ? = 1A8 in 3 ; I = 1100 in 4 

Z. = 67 in 3 ; A = 8,3 in 3 ; I. = A15 in 4 
4 4 4 



The retraining member sizes follow by symmetry, 
ties make up the matrix K, The subroutine PICK 
to obtain the properties of members one through 
routine COMSTRESS would be called to obtain the 
for members four through eight* 
p. Subr ou tine Stiff jaess 

i 

(1) Calculates new stiffness factors, k.. 

■ ij 

(2) Calculates the change in stiffness factors 
respect to the initial member choices. 



These proper- 
would be used 
three; sub- 
properties 



with 
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D. Subroutine Stiffness (Coat’d) 



(3) Determines AM and the new factitious maximum moment 



(4) 


Picks new 


nenbers 


based one 


This 


subroutine 


begins 


with the 


and 


calculates 


v 






t 




t 




k 17 = 


.310 


k 71 




k 17 = 


.184 


k 71 




A k ; 7 - 


+ .126 


Ak 71 



These results are used in equation (4) to obtain AM^ . 



AMj_ 7 = + 50.4 K f ; ll ^ - + 495 + 50.4 
Subroutine Pick now cones un with 



= 545 KIP FT 



z' = 153 in 3 ; I* = 1140 in 4 



These member properties are immediately inserted in the 
matrix K, This provides a constant updating of the data. 
The computer proceeds thru the structure doing the 
above for each member 
Member 2 



t » 



k 12 


.448 


k 21 


- .452 


k 12 = 


.450 


k 21 


= .450 


Ak 12 = 


-.002 


Ak 21 


= + .002 




= 1.1K’; 


" M 12 


= 635 KIP FT 



z' 2 - 204 in 3 ; I* - 1600 in 4 
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D. Subroutine St 


iffness (Cont 


’d) 


Member 3 






i 

k.. 


= .310 


t 

k = 0 


17 




71 


k,_ 


= .184 


k “ 0 


17 




71 


Ak 17 


- + .126 


£> 

H* 

It 

O 


f 




1 


> 

ti 


+ 50.4 KIP FT 


M ? = 545 KIP FT 


ii 

- r>. 

NJ 


3 1 

153 in' ; I = 


1140 in 4 


Member 4 








t 




1 


k.. 


= .117 


= 0 


13 


31 


k n * 


.184 


k =0 


13 


31 


» 




i 


Ak. _ 


- .067 


- 0 


13 




31 


f 




V 


AM = 


27 KIP FT 


M 3 = 73 KIP FT 


t 


1 ' 


L 


Z 3 = 


55 in ; = 310 in Subroutine COMSTRESS 


Member 5 


. i » 




t 


k !C = 


.121 


k c1 = 0 


15 


51 


Ak, _ = 


.184 


Ak c1 =0 


15 


51 


lk 15 ‘ 


-.063 


k 5 i = o 


I 




1 


a:i 5 = 


-25.2 KIP FT 


M 5 = 75 KIP FT 


V 


3 1 


2 ' 4 


Z 5 = 


55 in , A,. = 7 


,2 in ; Ij = 310 in 



Members 6 and 7 folio;*’ by symmetry . 
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The matrix K at this point becomes the matrix K . 
Subroutine Stiffness repeats the above process and produces matrix 

II 

K , The results of this procedure are summarized below 



Member 


' 3 

Z (in ) 


• 2 

A (in ) 


I (in 


1 


167 




1300 


2 


204 




1600 


4 


50 


6.9 


280 


5 


50 


6.9 


280 



The matrix K is now the matrix of actual member 

properties, A stress analysis of the structure made up of these 

members is performed; The results of this analysis are: 

M. = 595 KIP FT o, = 42 KSI - 1,07 o 
1 1 o 



M » 630 KIP FT 



M, = . 41 KIP FT 
4 



o _ = 40 KSI = o 
2 o 



a, = 36 KSI - , 9o 
4 o 



Stresses and moments for other members follow by symmetry. 
Subroutine Change now finds all the stress levels 

it 

satisfactory. The members of matrix K represent the optimum 
structure. The design spiral is complete, 

EXAMPLE 2 

This example was chosen to illustrate the procedure for 
a nonsynmetr ical system and one In which the field moment governs, 
at least initially. 
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The frame analyzed is shovn on page 46 . The optimization 
process proceeds as illustrated in detail in example one, 

A « I n i ti nl e r ^ P r o oe r t ies 

Z 2 = 240 in 3 ; I - 2000 in 4 

3 2 A 

All other members - Z = 40 in ; A = 6 in ; I = 210 in 
Jir- _ Stre s s Ana lysis 

M * l&S KIP FT 

M 2 =780 KIP FT 

M, =168 KIP FT 
4 

M = 84 KIP FT 

C_.__ Sub routine Change 



°1 * 


50 KSI 


1 

Z = 50 


3 ' 4 

in ; 1^ = 290 in 




°2 = 


3S KSI ; governing 

r 


z' = 240 


' in 3 ; I = 200 in 4 




°4 " 


50 KSI 




l 






Combined Stre 


O S | 


T I O 


1 


°5 = 


38.5 KSI 


Z. = 83 
4 


. 3 . n c . 2 

m ; A^ = 9,5 m ; 


*4 






Z 5 = 37 


3 ‘ ' 2 

in ; A = 5,6 in ; 


1 

X 5 


D, Subroutine 


Stiffness 








i 

Matrix K 










• i 

Member Z (in ) A 


• 2 

(in ) 


i’ a.. 4 , 




1 


56.1 


— 


320 




2 


203 


— 


1700 




4 


140 


13.5 


1050 




• 5 


23 


4.1 


100 





560 in 4 



= 190 in 
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The field moment governed in member 2 at this point. 
However once the second set of members were obtained, the end 
moments predominated. The results of the second pass is 
matrix K 



Member 

1 

2 

4 

5 



" 3 

Z (in ) 

56.1 

203 

190 

15 



" 2 
A (in ) 



16.5 

3 



i" (in 4 ) 
320 
1650 
1500 
60 



E. S tress Analys is 



M 1 = 530 KIP FT 

M 2 » 655 KIP FT 

M. = 114 KIP FT 
4 

m 5 = 0 

F. Sub rou tine^ Ch a nge 

o, = 38. 8 KSI = .97 

1 o 

o = 30 KSI = .75 

2 o 

o. = 48 KSI = 1.25 
4 o 



42 



These stress values require another pass. The trend is 
obvious, Members two, four and six are taking all the load. 
The second iteration with subroutine Stiffness produces the 
following matrices. 



i 

Matrix K - second pass 



! 3 

Z (in ) 

17 

203 

262 



Member 
1 
2 
A 

5 

If 

Matrix K - second pass 
Member 



" 3 

z (in ) 



1 

2 

A 

5 



192 

262 



(in 2 ) 



21 



A 



(in 2 ) 



21 



G, Stress Analysis 



I* (in 4 ) 
AO 
1650 
2000 



" A 

I Un) 



1550 

2000 



M 2 » 6A6 KIP FT 

= 6Ab KIP FT 

H, Subroutine Change 

o - 39.5 KSI = .98 a 

l o 

a, = 29.6 + 11.1 = A0. 7 KSI = 1.025 o 
A o 



- A3 



These stresses are within the specified limits, 

it 

The design spiral is terminated. The members of matrix K , 
second pass, form the optimum structure. 



12 


IU. 

© 

KIPS/FT 


AAA Ol 

3 

24 KIPS/FT 


AU 

4 

© AU 

L 

12 KIPS/FT 


L’s cqu 
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V 


* V 




.. ' 


f > 




\ 


/ 1 


s 
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© 




2 


© 




2 


© 
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K 






© 
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wv 


5 
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STRUCTURE FOR EXAMPLE B-l 



FIGURE B-l 
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STRUCTURE FOR EXAMPLE B-2 



FIGURE B"2 



APPENDIX c 



THE IDEAL SECTION 



The section moduli of available steel sections listed in 



the jal of^ St eel Con, struct i on do not vary in a continuous 
manner. The use of the least weight sections compounds the 
problem. Each section is useful for a wide range of section moduli. 



These facts make it very difficult, if one is using the 



sections listed in the AISC Manual, to pick out trends. Also the 
fact that there are no mathematical relationships existing between 
the important properties of the sections compounds the problems 
encountered when using actual sections. 



A mathematical section, i , e , one in which there are 



specific relationships existing between the section parameters, 
avoids these problems. The section developed for this thesis is 
outlined below. It conforms well to minimum weight sections listed 
in the AISC Manual for larger section moduli. At the lower end of 
the scale there is a greater divergence from actual sections. 



Let t^ = flange thickness 

t = web thickness 
w 

D = depth of section 
= flange width 

The relationships chosen for the section are; 



t 



f 



2 t 



w 



D = 50 t 



w 



w 






20 t 



w 
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The expressions for the 
of inertia are: 

A = .052 D 2 
Z = .032 D 3 
I - .0154 D 4 

A log-log plot of these 



area, section modulus and moment 



expressions is made on page 49* 
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BEAM DUI’TH (INCHES) 



FIGURE C-l 
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